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We prove that the gauge dependent gluon spin, gluon and quark orbital angular momenta oper- 
ators have gauge invariant expectation values on hadron states with definite momentum and polar- 
ization, therefore the conventional decomposition of nucleon spin into contributions from the spin 
and orbital angular momentum of quark and gluon is gauge independent. Similar conclusions apply 
to the gauge dependent quark momentum and kinetic energy operators, and accordingly nucleon 
momentum and mass structures. 
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qq' ' Gauge invariance is considered as one of the most fundamental principles in modern physics. All physical laws 

Q^^ , and observables are believed to be invariant under gauge transformations, and hence should be given gauge invariant 

On ' theoretical descriptions. Yet we sometimes do encounter gauge non-invariance. A case in point, which is well known in 

\ the contemporary extensive studies of nucleon spin structure, is that the gluon spin, gluon and quark orbital angular 

O I momenta operators are not separately gauge invariant. And in retrospect, the electron angular momentum operator 

^ whose expectation values are usually used to label the atomic energy levels, is not gauge invariant either. Then 
regarding the principle of gauge invariance, how should we, if not impossible, give definite physical meanings to these 

On operators, and especially to the labeling of atomic energy levels which we have talked abut for many decades? 

• One possible approach to solve this problem explored in the studies of nucleon spin structure is to find a particular 

\ gauge in which the gauge dependent operator can be related to experimental quantities. For example, in light-cone 

^ ' coordinates and light-cone gauge = the gluon spin can be related to the polarized gluon distribution measured 

\^ , in hard QCD processes In this paper, we will examine a more appealing possibility: since what are related to 

' experimental quantities are not the quantum operators themselves, but their matrix elements on physical states, we 

' will not worry about the gauge dependence of operators if they give gauge invariant matrix elements on the studied 
states. We find that this is indeed the case for the operators concerning hadron mass, momentum, and spin structures. 

QQ These operators, even though gauge dependent, have gauge invariant expectation values on hadron states with definite 

^\ momentum and polarization. 

To be specific, let's begin with the nucleon spin structure. The QCD rotation generator can be separated into four 

p I parts 

Ph! Jqcd = ^ y <f'xil)'^^i^ + j d^xifj'ix X 



-2 J d^xTr |i? X + 2 j (fixTr ^E^x x 



>< ■ =Sq+Lq+Sg+Lg, (1) 

' where and below any quantity is defined in the same way as in ref. and repeated indices are summed over. They 
have the obvious interpretations of quark spin, quark orbital angular momentum, gluon spin, and gluon orbital angular 
momentum respectively in the canonical formalism. However, there exists the well known difficulty that Lq, Sg, and 
Lg are not separately gauge invariant, therefore their physical significance seems obscure. Great efforts have been 
devoted to overcoming this difficulty. One example is what we have mentioned above for gluon spin, to study the 
experimental significance of these operators in a particular gauge. And on the other hand, one may naturally want 
to find gauge invariant operators for gluon spin, gluon and quark orbital angular momenta. Ji [Q and ourselves |^ 
have indeed obtained an explicitly gauge invariant separation of QCD angular momentum operator 

Jqcd = ^ J (fxip''f,ip + J d^xtp^x x -Dip 

+2 J d^xTr |x X [ex B^^ 
^Sq+L'g + J'g. (2) 

However, a further decomposition of J'g into gauge invariant gluon spin and orbital parts is lacking, and L'^ and 

do not obey the J x J = iJ algebra even when tp and A refer to bare fields, although Jg does in a pure gauge field 
theory. Therefore one should even hesitate a little to call them angular momentum operators. 
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Now go back to the gauge dependence of the decomposition in eq.(|T|). Does it really mean that we are unable 
to discuss gluon spin, gluon and quark orbital angular momenta without referring to a particular gauge? And 
consequently, only in a particular gauge does the decomposition of nucleon spin into corresponding contributions 
make sense? Or actually makes no sense at all since gauge invariance is committed so much priority in a gauge theory 
like QCD? A more realistic problem: have we made great negligence in the scores of years of talking about quark 
or electron angular momentum in quark models or atomic physics? And they actually also have no definite physical 
meanings? 

We find that the problem is, fortunately, not so serious. The point is that although the gluon spin, gluon and quark 
orbital angular momenta, or electron angular momentum operators are all gauge dependent, their contributions to 
nucleon or atomic spin, which are defined as their expectation values on nucleon or atomic states with definite 
momentum and polarization, are all gauge independent. We give below the proof. 



To obtain the matrix element like (piSi 



P2S2 ) (with \ps) a nucleon state with momentum p and polarization 



s), we consider the following three-point Green function: 



T 



[Ni{xi)LqNl{x2)} 



(3) 



where Ni{x) is the nucleon interpolating operator. 

As pi, P2 go to the nucleon mass shell, G^^l^S becomes 



Using Lorentz invariance to write 



'H\Jp\ + 



p\ + rn? 



le 



(27r)3^(0|A^;(0)bis) ipis 



P2S 



X {p2S 



nI,{Q) 0)(27r) 



-2i^p^ + ni^ 



V2 



(4) 



(0|iVz(0)|pis) = (2^)-3/2zu,(pis), 
P2s' A^/,(0) 0\ =(2^)-3/2z*u*(p2s'), 



where ui{ps) is the Dirac spinor with normalization '^iUi{ps)ui{ps') — S< 
{2tt)^^^^u'^{piSi){2tt)^^^^uii{p2S2) and summing over W , we get 



(5) 

and multiplying G^^^i ^ with 

^ I L q N ^ I 



IV 



p\ + vn? — ie 



Z\ (piSi 



P2S2 



-2i^Jp\ + m?- 



P2 + rn^ 



(6) 



The matrix element in eq.(|6|) is just what we want. The normalization factor \Z\ can be obtained from the 
two-point Green function 



NiN' 



d'^xid X2e 



ipi-Xl —ip2-X2 



q\t [Ni{xi)nI{x2)] 



-2i\Jp\ + m? 



p\ 



Z\'Y,MPisK'iPisK2n)HHPi+P2). 



(7) 



Eqs.(||, 0) show clearly that the gauge dependence of (piSi 



f-ypiSiP2S2 

NL,N 

to write 



P2S2 1 is equivalent to that of the ratio 



S2 iQPi^P2^^ . To explore the gauge dependence of G^^^^*^ and G^^^. , we use the method of functional integral 

'^Tz^nT = / rf'2;id4x2e-'P-^^e-*''--^^(2^)~3/2^;(p^,^)(2;r)-3/2^,,(p2S2) 
X J Y[dcl>Jx)e''Ni{xi)LgNl,{x2)B[m]DetT[cl)] 



(8) 
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(and likewise for G^^^'^ ^ ), where (j)n{x) denotes collectively the quark and gluon fields, B [/[</>]] is some functional of 
a general gauge-fixing functional fa [</>; x] , and JF is defined as 

J-ax,f3y[<P\ = ^A'3(y) '^=°' 

where (/)A(a;) is the result of after a gauge transformation with parameters A"(x). 

By saying that a matrix element or Green function is gauge independent, we mean that it is independent of the 
choice of the gauge-fixing functional and the functional B[f] up to an irrelevant field-independent constant 

factor. To demonstrate the gauge independence of G^i^^^^ and G^^^ , we first give a brief review of the standard 

proof of gauge independence for the Green function of gauge invariant operators O1O2 ■ ■ ■ Oi |^]: 

G01...0. = (0|T{Oi.--OJ|0) 

= ^lldMx)e''Ol■■■O^B[m]DctT [(!>]. (10) 

n,x 

Change the integration variable from cj) to (f)^, with uj"{x) an arbitrary set of gauge transformation parameters, and 
make use of the gauge invariance of the measure Yin x '^4>n{x), the action /, and the operators Oi • ■ • O^, we have 

Goi-o. = / n dcj)n{x)e"Oi - O^B [f[^^]] Det^F [<l>^] . (11) 

n,x 

It can be verified that Q 

DetJ^ [(f)^] = Dot J [0, Lu] ■ Detn [uj] , (12) 

where 

/SyH- g^isi^y-^ Ia=o (13) 

(with uj" the parameters of the associate transformation e^*'^(^)e~*'^'-^^ = e"*'^'-^^) is solely a functional of uj"{x), and 

Det^[0,^]=Det^:||^ (14) 

is just the Jacobian when transforming the integration variables from uj^iy) to fa['f'(^',x]. Therefore if we multiply 
both sides of eq.(|ll|) with the weight functional 



p[uj] = 1/Det7e M (15) 

and integrate over uj"{x), then by noting that Goi - Oi is independent of lo"{x), we get 

SI{n,xdK{x)e^'Ol■■■0^m 
(-^Oi - o, ppf J at \ \ ^ ' ^^^1 

J Ila,xd'^ i^)PM 

where 

C[</>] = / l[dco''{x)p[uj]B[f[cj)^]]DetT[(l>^] (17) 

a,x 

carries all the dependence of Goi - Oi on the gauge-fixing fmictional f[(f);x] and the functional B [/]. Now make use 
of eqs.(|l2|, we have 

cm = / l[duj''{x)BetJ[cl>,Lu]B[f[cl>^]] 

Cy.,X 

= J Y[dfa{x)B[f] = C. (18) 
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This is clearly independent of the choice of f[(j), x], which has been reduced to a mere variable of integration, and 
it depends on B [/] only through a (/)-independent constant C. Thus we proved that Goi - Oi is gauge independent as 
the operators Oi • • ■ Oi are gauge invariant. 

Now we come to G^t^*^^^ and G^j^^^^^ . Since nucleon is color singlet, the nucleon interpolating operator Ni{x) is 

SU(3) gauge invariant. Therefore according to the above proof, G^^^^ ^ is SU(3) gauge independent. However, Lq is 
SU(3) non-invariant. After an arbitrary gauge transformation 

'4>{x) u{x)4', 

A^ix) ^ Uix)A^U-\x) + -U{x)d^U-\x) (19) 



Lq induces an extra term 



5Lq = J d^xtl^l (^U-^x X ^dU^ V, 



(20) 



where U{x) — e — e ^'^ '^^^"^ . So if we also change the integration variable from (j) to as we did for Goi - Oi, 



QPisjp,s2 becomes: 

NLaN^ 



^nEnT = / d':^;id'2;2e-'f-^^e-^P--^^(2^)-3/2«;(piSi)(27r)-3/2y,,(p2S2) 
\{dcpn{x)e'' Ni{xi){Lq + 5Lq)Nl{x2)B [f[(t>u\]'DciF[^^] 

n.x 

^Gl,+G,z,. (21) 

Due to the extra term G^^ , the above proof of gauge invariance does not simply apply. However, if we can show 
that Ggj^ actually equals zero, then the above proof applies and hence G'^^^'jl^^^ is gauge independent. 
Comparing the expression of G^^ with eqs.(^, ^, ^) shows 

GsL,^ j d4a;id4x2e-'''--ie-'P-^^^(2^)-3/2u;(p^si)(2^)-3/2y;,(p2S2)(o|T{A^Ka;i)'5i;^^^^ 

2 ,^2_,v 1^1 (piSilSL^ P2S2) ^2 ^^2 ' (22) 



Pi + — le 

where 



P2 + — IS 



5L'q 



cPxi;^ (^Ux X ^dU-^^ V- (23) 



Without losing generality, we let the nucleon polarize along the third direction: JqcdIp^) ~ ^\P^)- We will 
demonstrate that SLq' can be expressed as a commutator with Jqqu, then eq.(22) shows clearly the result we aimed 
at: Gf-f vanishes if si = S2 and hence G^Jp^^^^^, is gauge independent. 

To obtain the commutator expression of 5L'^, we express the pure gauge form U(x)df^U^^{x) in terms uj{x) 



Define Ur = e^'^^ , V'r = f^rV', and note {x x d) = d^p, we get 



Ud^^U-^ = e-^^a^e*'^ = I I drer""^ {d^^Lo) e""^ . (24) 
3 



4 







= - f dT j d^'x {d^ {'^luoil^r) - Jd^ (l^t^Vr) } 

dT J d^xe'^^dj [xiiPlioi^r) 

- I dT [ d^xi[j^cD,^l^^r] , (25) 







where in the last step we have noticed that is a c- function commuting with Jqqj), and iplT"ipT is a scalar 
operator under spatial rotation, therefore 9^ ('iplT°'^pr) = —i [jQCDT''pT'T°''^r\ , which is analogous to the Heisenberg 
equation of motion dtO = i [Hqcd,0]. Another critical point in this observation is that, as can be seen from the 
explicitly gauge invariant expression in eq.(|2|), the total angular momentum operator Jqcd is invariant under gauge 
transformations, i.e., it can be equivalently expressed in terms of or ^jj^, AJ^ = U'^^A^jUt + ^U~^df^Ur- 

The first term in the final expression of eq.(^5|) is a surface term. Note that as \x\ — > oo, ip'l'4'T vanishes faster than 
ja?!""^, therefore this term vanishes after integration. Such surface terms will be encountered again and again later in 
this paper, and can all be shown to vanish by similar arguments as here, so we will drop them out without further 
explanations. 

Now all we are left with in eq.(|25|) is a commutator, as we aimed to show. Thus we have proved that both G^l?^^ 



and G^^^t are SU3) gauge independent. Then according to our above explanations following eq.(^, {pis 



P2S ) IS 



SU(3) gauge independent. 

In the case of U(l) gauge transformation, the pure gauge form Ud^U^^ is just id^uo, so some of the above algebras 
become simpler. However, since the proton carries electric charge e, Ni{xi)Nj,{x2) acquires an extra phase factor 
^-ieui{x^^ieuj{x2) after U(l) gauge transformation. But this phase factor makes no trouble: the only modification 
we need is to change the choice of weight-functional in eq. (||) to p[uj] = e*'='^(^i)e-*^'^("=^VDet7eH; and without 
modifying any of the above conclusions."'^ Therefore (pis Lq P2s) is also U(l) gauge independent, and we do not 



have to worry about the physical meaning of the labeling of atomic energy levels any more in spite of the gauge 
dependence of the electron angular momentum operator. 

We summarize our above analyses as follows: if after an arbitrary gauge transformation with parameters a;"(a;), 
the extra term SO^j acquired by some gauge dependent operator O can be expressed as commutators with the rotation 
and/or translation generators of the system, then O has gauge independent expectation value on the states with 
definite momentum and polarization. (Actually what appears in the calculations is SO'^, the result of applying an 
"anti-transformation" with parameters ~lu"{x) to SOui, as shown in eqs.(|2^, p3| ). But since applying this "anti- 
transformation" to (O + SO^) just gives the original O, 50'^ is just —50-^\ therefore if 50^ can be expressed as a 
commutator, so does 50'^.) In the remainder of this paper, we will show that all the operators concerning nucleon 
spin, momentum, and mass structures meet this requirement. 

First we look at the other two gauge dependent operators concerning spin structure: Lg and S^. Under the above 
gauge transformation in eq. (p^) , we have (from now on Ut is defined as e~''^" instead of e''^") 

5Ll^2 j d^xTr i^E,U-^d^ {uAiU-^ + V^.C/^i^ - E.d^Ai 

= 2 j^dT j d^xTr |i [El, Al] O^uj + i {^^El ~ Ur {d^E,) U~^) d^uj^ 



^This technique does not apply, however, if an operator is not SU(3) invariant. For example, after the gauge transformation 
in (|ie|), 'il^c,{xi)ilif3{x2) becomes Uaa'{xi)'ilJc,'{xi)-il>f3'{x2)Ufsifs{x2)- In such case we are unable to make a simple remedy by 
changing the choice of p[u}], and the above proof of gauge invariance fails. Of course, further investigations might show that 
such gauge dependence cancels in the ratio of three-point to two-point Green functions, leaving the studied matrix element 
still gauge independent. However, if this is not the case, then on a color non-singlet state, the expectation values of Lq, and 
even the expectation values of gauge invariant operators, will be gauge non-invanant; and we get a peculiar and surprising 
connection between gauge invariance and color confinement: gauge invariance requires color singlet. 
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9 Jo J 

^I^T J d^xTrS^id^ij) (^t [EJ,Al] + ^d,El^ - ^e'^^ EJ^ 



(26) 



Use diO — ~i [Pqcd^O] addition to d^pO — —i [Jqcd^O\^ vfiih Pqqjj the QCD translation operator (see 
eqs.(|3^, |3^) below), all the terms can be written as commutators: 



5Ll^2 I dr I d-x 



4cD,Tr[Lu[AlEl] + -cod,El 



-- f dr f d'xe'^' [Pl^cD:Tr{u;Ej)] 
9 Jo J 



(27) 



The calculation for 5^ is much simpler: 



^5*^ = ^^ dr J d^xe'^^Tr{{^^uJ)EJ} 

f^dr j d^xe'^^ [PhcD.Tr [loE])] . (28) 



2i 
9 Jo 



Thus according to our above theorem, we come to the conclusion that the quark's and gluon's spin and orbital 
contributions to nucleon spin, which are defined as the expectation values of S^, L^, Sg, and Lg on a nucleon state 
with definite momentum and polarization along the third direction, are all separately gauge independent, even though 
the operators L^, Sg, and Lg are not. 

It should also be noted that the expectation values of L^, Sg, and Lg are not separately Lorentz invariant. In 
the infinite momentum frame of the nucleon (or equivalently, in the light-cone formalism) and the light-cone gauge 
= 0, the first moment of the polarized gluon distribution measured in hard QCD processes can be related to the 
above local operator of gluon spin, which in the light-cone coordinates is accordingly Q 

S+ = J d^xMl'' 

= 2 I d^xTrl (e X aY + A± ■ bA , (29) 



where M^^'' is the angular momentum tensor of QCD and the subscript Sg denotes the gluon spin part. The 
expectation value of S'^ on the \ps) state is obviously also gauge invariant, since in our above proving of gauge 
invariance, the nucleon momentum is not constrained. 

The discussions for nucleon momentum and mass structures are similar. The QCD momentum operator can be 
written 

pQCD = J d^'x^P^^diP + 2 J d^xTr (^E.dAi) , (30) 
or in an explicitly gauge invariant form 

Pqcd ^ J d^x^j^D^j + 2 j d^xTr (e x (31) 

where Pq — J d^xip^ jd^p is the gauge dependent quark momentum operator. By the same techniques as above, 

it is easy to show that dPq can be expressed as a commutator with Pqcd, and hence the quark contribution to 
nucleon momentum is gauge independent. Since the total nucleon momentum is gauge independent, so does the gluon 
contribution. 

The QCD Hamiltonian operator can be written 



H, 



QCD 



d^x%l)^ • -D + mf]J V 

) , (32) 



d'xTr [E^ + B^ 
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where / (fixtp^a ■ ^Dip is the (gauge invariant) operator for quark kinetic and potential energy, and it was shown by 
Ji that it contributes about one third of the nucleon mass On the other hand, the quark kinetic energy operator 
Ek — J (Pxip'^d ■ jdtp is gauge non-invariant, but the extra term it acquires after a gauge transformation can also be 
expressed as a commutator 



5Ek = - y dT j d^x%l>\d ■ (^duj^ Vr 

= -^J^ dT J d^x [P^cD,^^a^^^] , (33) 

therefore the quark kinetic energy inside the nucleon is still gauge independent. 

In summary, we showed that gauge dependent operators are possible to have gauge invariant matrix elements 
on some special physical states, and thus can have definite physical meanings and experimental significance. This 
includes all the gauge dependent operators encountered in the studies of hadron spin, momentum, and mass structures, 
and also of atomic energy levels. These operators have gauge independent expectation values on hadron or atomic 
states with definite momentum and polarization. Therefore the conventional decomposition of nucleon spin into 
contributions from the spin and orbital angular momentum of quark and gluon is gauge independent, although some 
of the corresponding operators are not; and similar conclusions apply to hadron energy and momentum structures. 

We owe special thanks to T. Goldman, X. Ji, W. Lu, G. Sterman, and G.W. Wong for fruitful communications. 
This work is supported in part by the NSF (19675018), SEDC and SSTC of China. 

Note added: What we obtained in this paper is a non-perturbative proof of gauge invariance, while gauge invariance 
to each order of perturbative expansion does not necessarily follow, and needs separate studies. Some relevant 
perturbative calculations can be found in ref. M. 
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